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SUMMARY 

The detailed construction of the general solution of a second order non-homogenous linear operator- 
difference equation is presented. The wide applicability of such an equation as well as the usefulness 
of its resolutive formula is shown by studying some applications belonging to different mathematical 
contexts. 
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1 INTRODUCTION 

In this paper we report the explicit representation of the general solution of the second order 
non-homogenous linear operator-difference equation 

Yp+2 = £>oY p + £>iY p+ i + (frp+i, (1) 

where the unknown {YJ,} pe pj as well as the non-homogenous term {cp p } p€ ^ are sequences from a 
vectorial space V, and the coefficients £q,£i ; are linear noncommutative operators mapping V 
on itself, independent from the discrete variable p £ N. This equation encompasses interesting 
problems arising in very different scenarios. If, for instance, the reference space V is the complex 
Euclidean space C n , that is Y p and 4>p+i are ^-dimensional vectors, £jq and &i n x n complex 
matrices, then eq. ([T]) is the vectorial representation of a system of second-order linear non- 
homogenous difference equations. As another example, let's identify V as the vectorial space of 
all linear operators defined on a given Hilbert space. Now, the operators &o and H\ act upon 
operators and for this reason are called superoperators. The master equations appearing in the 
theory of open quantum systems provide examples of equations belonging to this class [5]. It is of 
relevance to emphasize from the very beginning that the ingredients Y p , (p p , Lq and £i of eq. ([1]) 
may be also interpreted as elements of an assigned algebra V. Let's consider, for example, V as the 
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noncommutative algebra of all square matrices of order n, that is M n [C]. Then, eq. ([I]) defines a 
second order non-homogenous linear matrix-difference equation, where Y p ,(f)p,£jQ and £1 belong to 
M n (C). We wish further emphasis that if V is the vectorial space of the smooth functions over an 
interval /, that is C°°(I), then eq. (p} represents a wide class of functional-difference equations [1], 
including difference-differential equations or integro-difference equations [2-4] . These few examples 
motivate the interest toward the search of techniques for solving the operator eq.(pQ), with &q,£ji 
noncommutative coefficients. 

In this paper we cope with such a problem and succeed in giving its explicit solution leaving 
unspecified the abstract "support space" wherein eq. ([1]) is formulated. This means that we do not 
choose from the very beginning the mathematical nature of its ingredients, rather we only require 
that all the symbols and operations appearing in eq. (pQ) are meaningful. Accordingly, "vectors" 
Y p may be added, this operation being commutative and, at the same time, may be acted upon 
by £o ° r &i ( hereafter called operators) transforming themselves into other "vectors" of V. The 
symbol Yq = simply denotes, as usual, the neutral element of the underlying space. Finally we put 
(£j a £jb)Y = £j a (Lj,Y) = £j a £jbY with a or b = 0, 1 and define addition between operators through 
linearity. 

The paper is organized as follows. 
The first section presents the solution of an arbitrary Cauchy problem associated with eq. ([1]). 
Some interesting consequences of such a result are derived in the subsequent section. The practical 
usefulness of our resolutive formula is shown in the third section where we solve some nontrivial 
functional-difference and integral-difference equations. Some concluding remarks are presented in 
the last section. 



2 EXPLICIT CONSTRUCTION OF THE RESOLUTIVE FOR- 
MULA OF EQ.® 



Thus, as for the linear differential equations, and independently from the noncommutative nature 
of £jo and £i, solving eq. (pQ) amounts at being able to construct the general integral of eq. ([2]) and 
to find out a particular solution of eq. ([I]). To this end, we start with the following theorem which 
extends a recently published result [7] concerning the exact resolution of the following Cauchy 
problem 



Let's begin by recalling that if {l^*} pg N and {YJ,} pe N are solutions of eq.([T]), then {Y p H } p ^ 
defined as Y p — Y* = Y p is a solution of the associated homogenous equation 



Yp+2 — £o^p + Cj\Y p+ i 



(2) 




(3) 



Theorem 1. The solution of the Cauchy problem 




(4) 



can be written as 




(5) 
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where the operators a p and (3 P have the following form 



( i p- 2 

2 



a,. 



E {^^ 2t) }^ 



t=0 



E 



if p>2 

if p=l 
if p = 



(6) 



1% 



E 

t=0 





(*) p (p- 1 - 2 *) 



z'/ p = 1 



(7) 



We recall that the mathematical symbol fL^}, in accordance with ref [7], denotes the sum 
of all possible distinct permutations of u factors Hq and v factors H>\, while 0, E : V — > V define 
the null and the identity operator in V, respectively. We omit the proof of this theorem since it is 
practically coincident with that given in ref [7]. Here instead we demonstrate the following 



Theorem 2. Eq. admits the particular solution 




if p>2 
if p = 0,l 



Proof. It is immediate to verify, by direct substitution, that the sequence given by eq. ( 
eq. (P) written for p = and p = 1. To this end, it is enough to exploit eqs. ([7]) and 
Y* = = fa and y 3 * = fofa + 1( j)2 = + 02- 

For a generic p > 2, introducing Y* in the right hand of eq. (pQ) yields 



(8) 

satisfies 
) getting 



p-1 



r=l 



r=l 



p-1 



r=l 



(9) 



Applying theorem (1) to the Cauchy problem expressed by eq. ([3]), we easily deduce that for p > 2 
and r = 1, 2, . . . ,p — 1 the following operator identity 



&O0p-r + £l/3p+l— r — Pp+2—n 
holds. Thus, the expression given by eq. ([9]) may be cast as follows 



(10) 



p— i p+i 

^/?p+2-r^r + /?p+2-(p)0p + /?p+2-(p+l)0p+l = ^ Pp+2-r4>r 



(11) 



r=l 



where we have exploited the identity = £iA based on eq. ([7]). Since the right hand of eq. (fTT|) 
coincides with K,* + 2 as given by eq. (jHJ), we may conclude that {Yp*} pe N, expressed by eq. (jSj), 
provides a particular solution of eq. ([!]) . □ 
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On the basis of theorem (1) and (2) we hence may state our main result, that is 
Theorem 3. The solution of the Cauchy problem 

Yp+2 = &oY p + &iY p+ i + 4> p +i 



Y = A, Y X = B ' (12) 

is 

p-i 

Y p = a p A + P p B + Y^ Pp-Ar (13) 

r=l 

where A and B are generic admissible initial conditions and a p and (3 P are defined by eqs. ^ and 
([7|) ; respectively. 

We emphasize that eq. ()13j) furnishes a recipe to solve explicitly, that is in terms of its ingredients 
£o,£i and {0 p+ i} pg N, the general Cauchy problem expressed by eq. (fT2|) , In the subsequent 
sections we will highlight that our result is effectively exploitable, providing indeed a useful approach 
to solve problems belonging to very different mathematical contexts. This circumstance adds a 
further robust motivation to investigate eq. (pQ) and its consequences. 

We conclude this section looking for the structural form assumed by eq. (|13p solely relaxing 
the noncommutativity between the two operator coefficients and L\. To this end, it is useful to 
recall the definition of the Chebyshev polynomials of the second kind li p (x),x G C [8] 

b/2] . v 

Indeed, taking into consideration that the number of all the different terms appearing in the operator 
symbol {£q"^£^} coincides with the binomial coefficient f u + v \ ^ w ^h m _ m i n ( Uj w ) as we u 

as assuming the existence of the operator (— £o) 5 , then the operators a p and f3 p for p > 2 may 
be cast as follows 

,2 

and 



-(-£ )%-2 ( -£j 1 (-£j )~2 ) (15) 



Pp = (--GoJ^Up-i Q-Ci(-^o) _ (16) 

where U p f§£a(— £0) means the operator value of the polynomial U p defined in accordance 

with eq. (|14p for x = ^£i(— Lq) . Thus, the solution of Cauchy problem (|12p may be rewritten, 
for p > 2, as 

Y p = (-L ) P -^Up^ (^i(-Zo)-^ B - (-£ )3U p -2 Q-dC-'Co)"^ ^4 + 

P- 1 /1 \ 

r=l ^ ' 

where Iq = ^- and Y\ = B are the prescribed initial conditions. 
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3 SOME CONSEQUENCES OF THE RESOLUTIVE FORMULA 



The mathematical literature offers several ways of solving linear second difference equations 
such as the matrix method or the generating function method. In the following we will heuristically 
generalized these methods to the operator case. The novelty of our method enables to deduce, by 
comparison with these approaches, some interesting consequent identities. Indeed, the second-order 
operator difference equation ([I]) may be traced back to the first-order vectorial representation 



where Y r 



Y 



P+i 



Ci 



E 



Yp+i — C\Y p + 



(18) 





<pp+l 



A 
B 



Successive iterations easily lead us to the formal solution 

v 



Y v = ClY^Y. C l r ^ 



(19) 



r=l 



On this basis, the solution of eq.([T|) may be written as [11] 



p 



Y P = P 1 C P 1 Y + P 1 VCf- r i?i$ 



r=l 



(20) 



where Pi = {E 0) and R\ 





E 



This solution is of practical use only if we are able to evaluate 



the general integer power of the companion matrix C\. Exploiting our procedure of writing the 
solution of eq. (JT]), the vector Y p may be expressed, accordingly with eq. (fT3|) . in terms of operator 
sequences a p and f3 p like 



p-i 

a p A + P p B + Pp-Ar 
r=l 
P 

a p+ iA + Pp+xB + 0p+i-i 

r=l 



\ 



J 



a,. 



CXp+l Pp+l 



A 
B 



+ 



I p-i 

/] Pp—rfir 
r=l 
P 

\ Yl Pp+l-r<t>i 
\ r=l 



(21) 



Confining ourselves to the homogenous version of eq. ((!]), that is putting p +i = into eqs. ([E 
and (|2ip . we get the formula for the p-th power of the companion matrix C\ as follows 



a,, 



ftp 



(22) 



E 

Lq L\ J \ a p+ \ Pp+x 

Another possible way of treating eq.([T]) is via the generating functions method [9,10]. We recall 
that, given the sequence {l^,} pg N, the associated generating function Y(s),s G C is defined as 



c, Yp = Y{s)=Y,Y pt 

p=0 



(23) 



under the assumption that the series converges when |s| < ^, for some positive number £. The 
advantage of this method consists in the systematical possibility of transforming a difference equa- 
tion in an algebraic one in the unknown Y(s). In order to apply such approach to the operator- 
difference equation given by ([1]), we stipulate that S[^i^p], ^i(SY p ) are both defined and that 
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S[^i^p] = &i(QY p ),i = 0, 1. Accordingly, heuristically, we transform both sides of eq. (p} getting 



Y(s)-A-Bs Y(s)-A 

—-^ — = L x ^^ + L Y(s) + $(s) 

s z s 

Thus, assuming the existence of (E — Lis — Lqs 2 ) -1 within the convergence disk, we have 

Y(s) = (E- L x s - Los 2 )" 1 A +(E- L x s - L s 2 y\B - £ Jl A)s + (E - Lis - L s 2 )~ 1 ^(s)s 2 (24) 
or equivalently 

Y(s) = {E- L x s - L s 2 )- l [E - L lS ]A + (E - Lis - Lqs^Bs + (E - Lis - L s 2 )- 1 ^>(s)s 2 (25) 
On the other hand, accordingly with eqs. (|23|) and (| 13 j) it holds that 

oo p— 1 

Y(s) = Y,( a P A + P P B + Y, P P -rA)s p (26) 
p=0 r=l 

Thus, one notes that imposing (fi p +i = and B = 0, respectively A = 0, we heuristically find the 
generating function of the operator sequences a p and f3 p in the closed form as 

oo 

Sa p = ^2a p s p := (E — Lis - LoS 2 )~ l [E — Lis] (27) 

p=0 

respectively 

oo 

SP p = ^2i3 pS P ^(E-Lis-Los 2 )-^ (28) 

p=0 

The particular case Lq = —E reproduces the generating functions of the Chebyshev polynomials 
of second kind. Extracting indeed for the sake of convenience the first two terms of the series, that 

oo oo 

is writing ^ a p s p = ao + ctis + Yl a p sP with the help of eq. (fT5j) and ao = E, ai = we get 

p=0 p=2 

oo „ 

(E-Lis + s 2 )- 1 [E-Lis]=E-J2^p-2[ -^]s p (29) 

p=2 ^ 

The eq. ([29]) easily determines the generating function of the sequence {U p [Li/2]} pe ^ in the form 

oo „ 

U P'^Y^ P = { ~ E ~ LlS + s2 ^ E -^ + s 2 -E + Lis] 

p=2 

oo „ 

^Y, u pIy ]sP ' = ( E ~ L i s + s2 y lE ( 3 °) 

p=0 

The novel results obtained in this paper exploiting our resolutive formula ( eqs. (f22l) . (|271) . (j28l) ) 
clearly evidence that our recipe to manage eq.(H|) successfully integrate with other resolutive meth- 
ods. Thus, we may claim that our resolutive formula do not possess a formal character only, since 
it helps to provide new interesting identities. 
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4 APPLICATIONS OF OUR RESOLUTIVE FORMULA 



An example of matrix- difference equation coped with our formula 

Let us consider the second order matrix-difference eqution 

Y p+2 = M Y P + MxY p+1 + % +1 (31) 

where Mo, Mi E M n (C) are noncommutative nilpotent matrices of index 2, that is Mf = 0, i = 
0, 1. Prescribing the initial conditions Yq = A, Y\ = B, then the solution of this equation is given by 
the eq. (|13H . The analysis of the matrix term {Mq U ^m[^} which appears in the composition of the 
matrix-operator a p and (3 P brings to light interesting peculiarities due to the specific nature of the 

coefficients Mq and M\. By definition, the term {Mq^M^} represents the sum of all possible terms 

of u factor Mq and v factors M\. Thus, it is quite simple to deduce that now the matrix-term of the 

form M^ 1 M" 2 Mq 3 M^ 4 ... is equal with zero, if Vi > 1, (V)i. Hence, we deduce that the operator 

{Mq^M^} survives solely when u = votv = u±1. Indeed, when u = v,u > 2, then in the 

sum survive only the terms [MoMi] [MoMi] . . . [MqMi] and [Mi Mo] [Mi Mo] . . . [Mi Mo] . 

v ' » , ' 

u— times u— times 

Further, the only nonvanishing matrix-terms {M (w ' ) M 1 (u+1) } are M x [M M 1 ][M Mi] . . . [M Mi], as 

V v ' 

u— times 

well as from {M^ ) M 1 (u_1) } the terms [M Mi][M Mi] . . . [M Mi] M , respectively. Exploiting the 

V v ' 

(it— 1)— times 

above results for the matrix term {M^ M^f 1 2t ^} we establish that 



{M^ 3 Wf 3 = [M Mi] . . . [M Mi] + [M X M Q ] . . . [MiM ], p = ?>k + 1 



k— times 



k— times 

{M y ~'M)~'} = Mi [M Mi] . . . [M Mi], p = 3k + 2 

V V ' 

k— times 



(32) 



{MpM[ P * 3) } = [MoMi] . . . [MqM x ] Mq, p = 3k, 



(fc— 1)— times 
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where k = 1, 2, . . . . Similarly, we have that 



{M^^M^^jMo = [M Q Mx] . . . [M Mi] M , p = 3k + 2 



k— times 



{M^ 3 ) m[* ) }M = Mi [M Mi] . . . [M Mi] M , p = 3k + 3 

v « ' 

(A;)— times 

{Mq^M^HMq = ([M Mi] . . . [M Afi] M )M = 0, p = 3k + 1 

" V ' 

(fe— 1)— times 



(33) 



Hence, we may write the solution into a closed form 



<5p-2 r £ ^2 1 (M Mi) 1 — J M + <5p-3 rp^siM^MoMi) 1 — J M 

3 4 3 J 3 4 3 J 



[<?g-i [E -i 1 (M Mi)[ £ 3 i ] +8 P -i r £ -i 1 (MiM ) [£ 3 i ] +5p-2 r£ _2 1 Mi(M M 1 ) [£ 3 2] + 5^ r P - 31 (MoMpt^Mrj-B 

3 >L 3 J 3 >L 3 J 3 4 3 J 3 4 3 J 



-1 



r=l 



p — r — 1 1 



+ 2^[^ P-r-i r p-r-i i (^(MoMi) 1 ^— J + (Mi M ) 



<5 p-r-2 rp-r-2 ] M 1 (M M 1 ) [ILJ 3- J - ] + 5 p-r-3 ! p-r-3 1 (M M] ) [Ht^I M ] $ r 
3 4 3 J 3 4 3 J 



(34) 



j4n example of functional- difference equation coped with our method 



The three-term recurrence relation 



fp+2(t) = -f P (t - T ) + fp+l(t + n) 



(35) 



with the initial conditions A = /o(i) and -B = f\{t) is an example of functional difference equation, 
traceable back to eq. (JTJ). It is indeed well-known that if f(t) is a function of class C°°, then the 
translation of its independent variable from t to t + r can be represented as the effect on the same 

oo 

function of the operator exp[r4r] = ^ pbs]* - This operator appears in a natural way when 



fc=o 



one studies problems characterized by translational invariance in a physical context [12]. Thus, 
by putting £j = (— l)* +1 exp[(— l) l+1 Tj^], i = 0, 1 the commutativity property of the two operator 
coefficients £q and L\ allows us to write down the solution of eq. (|35p as follows 



f p (t) = exp[-(^-i)T ^]Up_i 



1 u . t , d, 

2 exp[(ri + y } ^ 



ex P[-(|) T 0^]^p-2 



1 u . To, d 
r x P [(r 1 + -)-] 



hit) - 
fo(t) 



(36) 
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Exploiting eq. (|14p we may write down that 

b-1/2] 



f p (t) = 

k=0 
[p-2/2] 

E (-D 

k=0 



p — 1 — k 

k 



/i (i + (p - 1 - 2fc)n - A;r ) 



p-2-k 
k 



f (t + (p - 2 - 2fc)n — (fc + l)r ) 



Imposing, for example, the following initial conditions /o(t) = e * and = e* we get 



/ p (t) = exp[-(^— -^)r ]ltp_i 



i ex P [(n + y )] 

-exp[-(n + y)] 



(37) 
(38) 



(39) 
(40) 



An example of integro- difference equation coped with our method 

Consider the difference-differential equation 

fp +2 (t) = f3f p+1 (t)+af p (t), a, (3 gR, p = 0,l,... 



(41) 



where is a C°°{I) function with /o(t), and {/ p (0), p = 0, 1, . . . } prescribed functions. 
The above equation may be rewritten in the equivalent form 



fp+z(t) = £i/p+i(*) + Lofptf) + /p+2(0) 



(42) 



where £o = £i = f3L and £(•) = J -dr. Eq. (j4"2l) is a particular case of eq.([I]). 

o 

The explicit solution of this equation requires the knowledge of the operator terms a p and f3 p . One 
remarks that f3 p is the sum of + 1 operator terms of the form {XLq/L^ -1-2 *}. Because £o = 

and ZLi = (3L then, for a finite p, holds 



{H k Q L p - 1 - 2k } = a k ^~ 1 ~ 2k ( P- l -' { 



LP- 



-l-k 



(43) 



Therefore, by direct substitution into eq. (|T|) it follows that 



Similarly, we have that 



k=0 



f P-2 



2/r / p — I — k \ „ p _l_/. 



(44) 



fc=0 ^ 



p — 2 — k \ ^p-i-fc 



(45) 
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The solution of the corresponding homogenous equation in accordance with the prescribed initial 
conditions is then 



f^® = E a fe+ V~ 2 - 2fc ( P " I " k ) (kit)) + 

^ a fe /3 p " 1_2fc ( P ~l~ k ) 



(46) 



fc=o 

Exploiting our central theorem (2), we may claim that 



p-1 p-2 

f; = £ /3 p _ m / m+ i(0) = ^ /3 P _ m / m+ i(0) + / p (0) (47) 

m=l m=l 

is the particular solution of the nonhomogenous equation for which /o = /i = 0, (V)t. Equivalently, 
we have that 

/p = E E a k /y- m - 1 - 2k ( P ~ m ~ 1 ~ k ]^- m - 1 - fc (/ m+1 (0)) + / p (0) (48) 

m=l fc=0 V / 

But, as shown in the Appendix, we may prove that 

t t n t n -i t:j ti t 

Z n (f(t)) = I dt n J dt n -! J dt n _ 2 ■■■ J dt 2 J f(h)dh = J(t- r) n ~ l f(r)dr, (49) 



so that we may write down that 

[2=1] t 

4 H \t) = Y: a fc+1 /3- 2 " 2fc ( P ~ l ~ k ) {p _l_ k)[ j (t - rY-^h{r)dr + 

£ ( ? ~ 1 _ " ) jpzhj^, ) (' - rr~AM* (50) 



fc=0 

and 



[ p-m-l 
P-2 [ 2 



/ P *=E E ^" m " 1_2fe ( P V * ) (p _ m 1 2 _ /(* - rr- m - 2 - k f m+1 (0)dr + / p (0) 
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Hence, the general solution of the proposed integral-difference equation is 



f^l t 



f p (t) = £ a ^f3^ ("'I'") {p _\_ k)[ J (t ~ Tf- 2 ~ k MT)dT + 

[2=1] * 

£ - k ^- 2k ("'I'") (p-2-fe), / (* - -r 2 -Vi(r)dr + (51) 



fc=0 

[ p — m — 1 1 

E E ( * - " - - ) _ m _ 2 1 _ m _ 1 _ + /,(0)(52) 

5 CONCLUSIVE REMARKS 

The novel and mean theoretical result of this paper is expressed by theorem (2) with which 
we demonstrate that eq. ([8]) provides a particular solution of eq. ([T]). This result together with 
theorem (1) completes the resolution of this equation enabling us to write down formula (|13p for 
its general solution. The operator character of eq. ([I]) and, as a consequence, the presence of 
generally noncommuting coefficients is the key to understand why such an equation may represent 
the canonical form of equations seemingly not related each other. The consequences of eq. f|13|) and 
the applications reported in this paper, besides being interesting in their own, demonstrate indeed 
both the wide applicability of eq. ([1]) as well as the practical usefulness of its resolutive formula. 

A APPENDIX 



For the sake of completeness we here report a proof of the well-known follwing identity 

t t n t n —l t2 t 

T)^f(T)<lT. (53) 



where f(t) is a C°°-function. The mathematical induction procedure will be exploited. 
For n = 1 the above formula becomes the identity 

t t 



t t„ l-n—l £3 12 C 

J dt n J dt n -i j dt n - 2 ---J dt 2 J f(ti)dt! = T^Tnj / (* 



J f(h)dh = J f{r)dr (54) 


Let's suppose that the formula [53] holds for any r < n and we prove it validity for n + 1. To proceed, 
it is convenient to define 

t t r t r — i £3 t2 

F r(t) = J dt r J dt r -i j dt r ^ 2 ---J dt 2 J f(h)dh (55) 


writing down eq.([53|) as follows 

F n (t) = ¥ ^J(t-TT-lf(T)dT 







t 
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It is obviously that by definition 

t 

F r+ x(t) = J F r (t r+1 )dt r+1 , r > 1 


By the induction hypothesis we have that 



(56) 



K+l - F n(t) 



(t_ T )«-V( T )dr 



(n- 1)! 

u 

Our problem becomes the resolution of the following Cauchy problem 



(57) 



(58) 



From the well-known Leibnitz identity 

b(x) b(x) 



^ j f(x,y)dy 



d f(x,y) 
dx 



F n+1 {0) = 



dy + b'(x)f{x, b(x)) - a'(x)f(x, a(x)) 



(59) 



a(x) 



a(x) 



we have that 

t 

1 



(n- 1)! 



(t- r )»-V(r)dr 



1 



o 



n(n — 1)! 



d[(t-rrf(r)} dr _ d 
dt dt 



o 



- I(t-rrf(r)dr 



(60) 



Hence, the differential equation for F n+ i(t) may be rewritten as 

t 



K+i(t) 



d_ 

dt 



1 



(t-rrf(T)dr 



o 



(61) 



such that the Cauchy problem (|58p has the solution 



F n+l (t) = - j{t-r) n f(r)dT 



(62) 
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